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Plasmons in spin-polarized graphene: A way to measure spin polarization 
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We study the collective charge excitations (plasmons) in spin-polarized graphene, and derive ex¬ 
plicit expressions for their dispersion in the undamped regime. From this, we are able to calculate 
the critical wave vector beyond which the plasmon enters the electron-hole continuum, its qual¬ 
ity factor decreasing sharply. We find that the value of the critical wave vector is strongly spin 
polarization-dependent, in a way that has no analog in ordinary two-dimensional electron gases. 

The origin of this effect is in the coupling between the plasmon and the inter-band electron-hole 
pairs of the minority spin carriers. We show that the effect is robust with respect to the inclusion of 
disorder and we suggest that it can be exploited to experimentally determine the spin polarization 
of graphene. 

PACS numbers: 73.20.Mf, 71.10.Ca, 81.05.ue, 71.45.Gm 


INTRODUCTION 

Collective oscillations of the electronic charge (also 
known as plasmons), after being for many years a topic 
of fundamental research in solid state physics [1-4], have 
recently emerged as a vector of information in a new gen¬ 
eration of electronic devices, compressing the energy of 
electromagnetic field to the nanometer scale [5-7]. At 
the same time, a major effort is underway to couple 
the spin of electrons to electric fields with the objec¬ 
tive to perform memory and logical operations. In this 
context, the interaction of plasmons with spins has re¬ 
ceived relatively little attention [8-10]. It has gener¬ 
ally been assumed that spin polarization, occurring ei¬ 
ther spontaneously, or by the application of a magnetic 
field, or by spin injection under non equilibrium condi¬ 
tions, would have only a minor effect on the propaga¬ 
tion of the charge mode. In a recent paper, we have 
pointed out the possibility that a new collective mode - 
a spin mode - might appear in a strongly spin-polarized 
Galilean-invariant two-dimensional electron gas (2DEG) 
within the pseudo-gap that separates the high-frequency 
single-particle excitations of the majority-spin electrons 
from the low-frequency single-particle excitations of the 
minority-spin component [10]. While this is a qualita¬ 
tively new feature of spin-polarized electron liquids, it 
does not affect the main characteristics of the charge plas¬ 
mon. 

The situation is profoundly different and far more in¬ 
teresting in the two-dimensional spin-polarized electron 
gas in graphene [11-15], which is the subject of this pa¬ 
per. At the root of the difference is the existence of 
low-energy inter-band excitations in which an electron 
from the lower cone of the massless Dirac fermion band 
structure is promoted to the upper cone without spin 
flip. These inter-band excitations are absent in the ordi¬ 
nary Galilean-invariant 2DEG. As the spin polarization 
of the electron gas increases, the energy of inter-band ex¬ 
citations of the minority-spin component decreases and 


eventually it becomes so small that it matches, at a crit¬ 
ical wave vector qc, the energy of the charged plasmon 
at the same wave vector. When this happens, the plas¬ 
mon ceases to be a well-defined excitation, since there is 
nothing to prevent it from decaying into the inter-band 
electron-hole (e-h) pairs of the minority-spin component 
(Landau damping). Experimentally, the onset of this de¬ 
cay process should show up as a sharp drop in the value of 
the quality factor, that is to say, the number of plasma 
oscillations that occur during an average plasmon life¬ 
time. This sharp onset, being strongly spin-dependent, 
offers a new and unexpected way to determine the spin 
polarization of electrons in graphene from measurements 
of the plasmon lifetime. 

This paper is organized as follows. In Sec. we re¬ 
view the RPA theory of plasmons in clean spin polarized 
graphene. We calculate the undamped plasmon disper¬ 
sion as a function of spin polarization and show that, at 
variance with the Galilean-invariant 2DEG, there can be 
no spin plasmon for any spin polarization, because the 
Fermi velocities of up- and down-spin electrons are iden¬ 
tical. We go on to compute the onset of Landau damping 
of the charge plasmon and the spin polarization depen¬ 
dence of the quality factor of this mode. In Sec. , we show 
that our results are robust against the inclusion of dis¬ 
order. Specifically, we show that short-range impurities 
in a high-mobility sample introduce a lower wave vector 
q* below which the charge plasmon becomes unstable to 
diffusion. Above this wave vector the plasmon disper¬ 
sion calculated in Sec. is essentially unchanged, and so 
is the sharp onset of Landau damping at q^. Section 
summarizes our results and conclusions. 


RPA THEORY OF PLASMONS IN GRAPHENE 

We consider doped monolayer graphene at zero- 
temperature, within the regime of applicability of the low 
energy linear-band dispersion relation [16], with a finite 


2 


spin polarization. Spin polarization in graphene can be 
induced through various means, such as non-equilibrium 
spin injection [17, 18], or optical excitation [19], or prox¬ 
imity effect [20]. However, for the sake of modeling spin- 
polarized graphene in the simplest possible way, we as¬ 
sume that the polarization is induced by an in-plane Zee- 
man field B, which has negligible orbital effects. For a 
given valley, the energy dispersion of the spin-resolved 
bands is given by 

ec<T(fc) = C + crghfJ-BB , (1) 

where ^ = -|- (—) indicates the conduction (valence) 
band, cr = -|- (—) indicates f-spin (4,-spin) electrons, vp is 
the Fermi velocity in graphene, which is typically 10®m/s 
— independent of the applied Zeeman field, and gb is the 
Lande g-factor, while is the Bohr magneton. The total 
density of states in graphene is given by (?s(;«jej/(27r?;|), 
with Qs = 2 and gy = 2 being the spin and valley de¬ 
generacies, respectively. The Fermi momentum in unpo¬ 
larized graphene is kp = yf^Knfg^y^ and the associated 
Fermi energy is £p = hvpkp. The spin polarization is de¬ 
fined as P = (n-f-n^)/(n-f+n^) « -sgn{B){g\,iJLhB/ep)'^ 
for small magnetic fields, where ria is the spin polarized 
carrier density in graphene. The spin-resolved Fermi 
wave vectors kpu can be expressed in terms of P as 
kpa- = kp\/l + aP. In the rest of the paper we will 
consider P > 1 without loss of generality. The Wigner- 
Seitz parameter which indicates the relative strength 
of Coulomb interactions, is a constant and it is given by 
Ts = /nhvp, where k is the static dielectric constant of 
the substrate hosting graphene monolayer. 

Following Refs. [2, 10] the collective density excitations 
within RPA in a spin-polarized medium are given by the 
zeros of the complex longitudinal dielectric function 
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FIG. 1: (a) Real and imaginary parts of the RPA dielectric 

function e{q,uj) vs ui for a 2DEG (parabolic dispersion) sys¬ 
tem. The zero of the dielectric function marked by the small 
circle corresponds to the spin plasmon mode in 2DEG [10]. 
(b) Real and imaginary parts e(q, uj) vs lj for graphene. Note 
that there is no spin plasmon mode in graphene. Other pa¬ 
rameters are chosen to be P = 0.5, q/kp — 0.1, and = 2.5 
for both panels. The thin vertical lines in both panels mark 
the location of the known plasmon mode in the long wave¬ 
length limit, i.e., in graphene utpi/ep = \j2raqlkp and in 
2DEG oipi/ep = 


given by 


7osc(g) 


STO[e(g,a;)] 
d^e [e{q,uj)]/duj 


( 3 ) 


In an actual experimental scenario, the physical observ¬ 
ability of a damped collective mode depends on the 
sharpness of its resonance peak, which is captured by 
the quality factor of the collective mode, defined as 
Q ~ ^osc{g')/^osc{g}• 

The spin-polarized Lindhard function for graphene [12] 
is best expressed in terms of the dimensionless variables. 


V±a- 



UJ 

vpq ’ 


( 4 ) 


e{q,uj) = 1 - Vg[xot(9,w)-bxot(g,w)] 

= 0 , ( 2 ) 

where xoa- are the spin-resolved response functions of 
the non-interacting spin-polarized medium, and Vq de¬ 
notes the Fourier transform of the screened or unscreened 
Coulomb potential. The oscillation frequency of vari¬ 
ous collective modes Wosc (typically charge and spin plas- 
mons) and the corresponding damping rate (or inverse 
lifetime) 7osc, is obtained by solving for the complex roots 
of Eq. (2): uj = Wosc(g) - *7osc(g), for a given q. Note 
that for the stability of the collective mode it is essen¬ 
tial that 7osc > 0. The collective mode typically decays 
by creating single electron-hole pairs (Landau damping), 
which occurs if the frequency of the collective mode lies 
in the e-h continuum, which in graphene can arise ei¬ 
ther from intra-band or inter-band e-h excitations. If the 
damping rate 7osc is small, it can be estimated by do¬ 
ing a Laurent-Taylor expansion of the dielectric function 
around ujosc [10]. For a frequency independent Vq, it is 


and at T = 0, for complex frequencies it is explicitly 
given by [15, 21] 


^ - G+{l'+yy)+G-{l'-a) , (5) 

Nq kp l J 

where 


F{q,uj) = 


Vpq" 


ikp^vlq^ - 


( 6 ) 


and 


G±[z) = z\/l — ± 7cosh ^{z) . (7) 


In the long wavelength limit, g —>■ 0 with uj > vpq 
fixed, which is important for optical spectroscopy and the 
plasma oscillations, the spin-polarized Lindhard function 
can be approximated as 


X0<T 

Nq 


{l + aPfG 





— v^q^ 


( 8 ) 
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As a check of Eq. (8) we note that a similar approxi¬ 
mate expression for the Lindhard function (without spin 
polarization) has also been derived in the context of topo¬ 
logical insulators, [22] which have the same quasiparticle 
dispersion as graphene. 

Having calculated the Lindhard function we now focus 
on the zeros of the dielectric function, i.e., the collec¬ 
tive modes. In Fig. 1, we plot the dielectric function of 
graphene [panel (b)] and compare it with the same for 
2DEG [panel (a)]. Evidently, there are two stable collec¬ 
tive modes in 2DEG, corresponding to charge and spin 
density excitations [10], while in graphene there is only 
one stable collective mode and it corresponds to charge 
plasmons. Physically this is a direct consequence of the 
fact that in a 2DEG with parabolic dispersion relation, 
the Fermi velocity of the majority spin carriers can be sig¬ 
nificantly larger than the Fermi velocity of the minority 
spin carriers. This opens up a region of reduced spectral 
density (pseudogap), which is clearly visible in Fig. la) 
between the peaks associated with e-h excitations of the 
majority and minority spin carriers: it is in this win¬ 
dow that the spin plasmon lives. At variance with this, 
the Fermi velocities of graphene are the same for both 
spin species and, as a consequence, there is no spectral 
window in which the spin plasmon can find a home (see 
Fig. 1). 

We now proceed to calculate the dispersion of the 
charge plasmon in spin-polarized graphene [11-14]. Sub¬ 
stituting Eq. (8) in Eq. (2) leads to the following approx¬ 
imate plasmon dispersion. 


Wpi(g) = v-pq 


1 + VgNodjP) 

^1 + 2VgNoS{P) ’ 


( 9 ) 


where we have defined S{P) = P + Vl — P. Note 
that S{P) is a monotonically decreasing function of P, 
varying from (5(0) = 2 to (5(1) = For the usual 

case of long range unscreened Coulomb repulsion in 2D, 
Vq = 2TTe^lKq, where k is the dielectric constant of the 
substrate, we have NoVq = 2rskp/q (with the factor of 
two originating from the valley degeneracy) and Eq. (9) 
can be rewritten as. 


Wpi((7) = vpy/q 


25{P)rskp + q 
yjA5{P)rskY + q 


( 10 ) 


The plasmon mode in spin polarized graphene enters the 
e-h continuum of the inter-band transitions of the mi¬ 
nority spin species when it crosses the region specified 
by w -I- vpq > 2A:FnF(l — PY^'^^ and the critical wave 
vector is obtained to be 

, A(P) + v/A(P)2 + 16fc2F,J(P)(l-P)3/2 

=- i7r=P -• 

( 11 ) 

where A{P) = fcF[(l — P) — 4rs(5(P)\/l ~ P — ''’s'^(-P)^]- 
In the limiting case of P —>■ 0, i.e., unpolarized graphene. 
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FIG. 2: The color plot of the loss function, —Am[l/e(q, ta)], 
along with the numerically evaluated plasmon dispersion 
(solid black line), and the approximate analytical dispersion 
(green circles) as per Eq. (10). Panels (a), (b), (c), and (d) 
correspond to P = 0, 0.25, 0.5, and 0.75 respectively. Here 
Ta = 0.5 in all the panels. 


Eq. (11) reduces to 


9c(0) 


kp 

Y 


1 - 8rs 


4rg -I- (1 -b 2rs)^/YYYrsYYA 


( 12 ) 

As a consistency check we note that Eq. (10) can be ex¬ 
panded around q —>■ 0 to obtain 


^pi(g) 

kpvp 




1/2 


-b 


8./YS{P) 



(13) 

Here the first term in Eq. (13) gives the well known 
behavior of the plasmon dispersion in 2D, and the second 
term, which also depends on the spin polarization, gives 
the additional correction which was reported earlier in 
the context of plasmons in unpolarized (P = 0, S{P) = 
2) intercalated graphite [11]. 

In Fig. 2 we plot the loss function and the plasmon dis¬ 
persion, both numerically calculated (solid black curve) 
and evaluated from the analytical expression (10) (green 
circles) for P = 0, 0.25, 0.5, and 0.75. Observe how 
increasing the degree of spin polarization leads to a de¬ 
creasing value of (/c, the critical wave vector for which 
the plasmon mode enters the inter-band continuum of 
the minority spin species, i.e., ujp\{qc) = 2?ivpkp^ — vpqc 
for P > 0. This is also reflected in the inverse of the 
quality factor as a function of q, which should be zero 
for q < qc since the mode is completely undamped in the 
present approximation, and will become finite as soon as 
the plasmon mode enters the inter-band e-h continuum 
of the minority spin carriers. This is clearly shown in 
Fig. 3 which also highlights the fact that the critical qc 
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FIG. 3: (a) vs g for various values of P in the clean 
limit. A sharp rise occurs when the charge plasmon enters 
the inter-band e-h continuum of the minority spin species and 
is approximately given by Eq. (11) and marked by thin verti¬ 
cal lines. A second cusp occurs when the plasmon enters the 
inter-band e-h continuum of the majority spin species, whose 
approximate location is specified by Eq. (11) with P —>■ —P, 
and becomes visible for large values of P. (b) shows the crit¬ 
ical qc beyond which the plasmon mode enters the inter-band 
e-h continuum. The P dependence of the critical qc arises al¬ 
most entirely from the shift in the boundary of the inter-band 
continuum: this is shown by the curve labeled “interband”, 
in which the polarization dependence of the plasmon disper¬ 
sion has been neglected. Further the maximum deviation of 
Eq. (11) from the exactly calculated critical wave vector oc¬ 
curs for P —>■ 0, and is less that 6.7%. We have taken Vs = 0.5 
in both the panels. 


decreases with increasing spin polarization. This jump in 
the quality factor of the plasmon mode aX q = Qc can be 
used to estimate Qc experimentally, and comparing this to 
either Eq. (11) or the exact numerical value of Qc, should 
give a good estimate of the degree of spin polarization P. 


EFFECT OF DIFFUSION ON THE PLASMON 
DISPERSION 


In this section we study the effect of weak disorder on 
the charge plasmon mode [23] in spin-polarized graphene. 
For sufficiently dilute impurities, the collective modes of 
the disordered system are still given by Eq. (2), pro¬ 
vided that the bare Lindhard function is replaced by the 
disorder-averaged response function. A suitable expres¬ 
sion for the latter was given by Mermin [24], based on a 
relaxation time approximation. It is given by 


noa(g,w) 


{uj + iT ^) xoa(g,^ + VT) 

UJ + iT-^X0cr{q,i^ + i/T)/xOa{q,0) 


(14) 


where xoaiq,!-^) is the Lindhard response function of the 
clean system, and r is the elastic lifetime of the momen¬ 
tum eigenstates in the presence of static disorder. Note 
that Eq. (14) is actually equivalent to the sum of impu¬ 
rity ladder diagrams (diffusons) in the diffusive regime 
[23] , and it also has the correct high-frequency behaviour 
(the collisionless regime). 

An approximate disorder-averaged Lindhard function 
can now be obtained by substituting Eq. (8) in Eq. (14), 




FIG. 4: (a) and (b) display Wpi, and Q~^, respectively, as 
functions of q for different disorder strengths, for P = 0.5. 
The shaded region around tUpi indicates the decay rate 7 pi, 
and the plasmon curves for different disorder strength are 
shifted vertically (by 0.4) for better visibility, (c) and (d) plot 
Q~^ for different polarization strengths, with 1/r = 0.02, for 
large- and small q, respectively. In (c) the cusps in the in¬ 
verse quality factor arise for the critical wave vectors where 
the plasmon mode enters the inter-band e-h continuum of the 
minority (first cusp - marked by vertical lines), and then the 
majority spin carriers (second cusp), (d) highlights the ex¬ 
istence of the P-dependent critical wave vector of Eq. (18), 
below which the plasmon modes are unstable to diffusion. No¬ 
tice that we plot Q rather than Q~^ and the different curves 
are shifted vertically (by 0.2) for better visibility. The sym¬ 
bols indicate the exact numerical result, which are in excellent 
agreement with the approximate analytical result of Eq. (19) 
indicated by the solid black line. In all panels, we have taken 
r, = 0.5. 


yielding 

« Vl + crP (-1 + 77-rwwT — ^2 211/2 -^ 

Nq \ [{uj-\- i/tY — v^q^y/^ — i/t 

(15) 

Note that, except for the spin dependent pre factors, 
Eq. (15) has a functional form similar to that of the dis¬ 
ordered response function of a 2DEG with parabolic dis¬ 
persion relation, evaluated by means of a perturbation 
theory in Ref. [23]. Now even the low energy plasmon 
dispersion will have an imaginary component and can be 
obtained by substituting Eq. (15) in Eq. (2). It is explic¬ 
itly given by 



1 + VqNo5 
1 + 2VqNo5 


^Jvlq^[l + 2VqNe)S\ - l/r2 


f/r , 
(16) 


where 5{P) is defined below Eq. (9). For unscreened 
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Coulomb repulsion Vq = 2'Ke^ jq, Eq. (16) reduces to 

\Jvlq[q + ^r^kpS] - -ijr . 

(17) 

Note that Eq. (17) reduces to Eq. (10) in the limiting 
case of vanishing disorder l/r —0. 

Equation (17) immediately implies that in spin- 
polarized graphene with dilute disorder, collective den¬ 
sity excitations can exist only for wave-vectors greater 
than a critical wave-vector (g*) which is analytically 
given by 

ql = -2r,kFS{P) + ^y{2rskp6{PW + ■ (18) 

Equation (18) is consistent with the physical intuition 
that coherent collective excitations cannot exist for 
length scales longer than the mean free path, i.e., long- 
wavelength plasmons “diffuse” in the disordered medium. 
The approximate value of the quality factor (for q > q^) 
can be obtained from Eq. (17), and is given by 

Q = T^vlq{q + ArskY5{P)) - l/r^ . (19) 

In panels b) and c) of Fig. 4 we plot the inverse qual¬ 
ity factor Q~^ vs q for various strengths of disorder (b) 
and degrees of polarization (c). The overall behavior of 
Q~^ is a U-shaped curve - the quality factor decreas¬ 
ing sharply both at large wave vectors (where Landau 
damping from inter-band excitations begins) and at low 
wave vector, where diffusion damping takes over. In the 
region of very small wave vectors, a plot of the quality 
factor itself (panel d) more conveniently illustrates the 
strong damping of the plasmon below the lower critical 
wave vector q*. The dependence of the lower critical 
wave vector q* on spin polarization could also be used, 
in principle, for an experimental estimation of P, but 
this approach would suffer from much larger uncertain¬ 
ties than the one based on the determination of the upper 
critical wave vector q^. 

CONCLUSION 

In summary, we have studied the collective density ex¬ 
citations of spin-polarized graphene, and derived explicit 
expressions for their dispersion and the critical wave vec¬ 
tor for the onset on Landau damping, considering the 
clean case as well as the weakly disordered case. Our an¬ 
alytical results are in excellent agreement with the exact 
numerical results. We have found that, unlike the case of 
2DEG, both the plasmon dispersion, and the upper crit¬ 
ical wave-vector beyond which the plasmon mode decays 
into the e-h continuum, are strongly spin polarization de¬ 
pendent. Further we have shown that this effect is robust 
with respect to the inclusion of disorder. The dependence 
of the upper critical wave-vector on spin polarization will 


manifest as a cusp in the measured quality factor, which 
we propose can be used to experimentally determine the 
spin polarization of graphene. 
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